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THE EFFECT OE ACCELERATING A HYPOTHETICAL AIRCRAFT 
THROUGH TEE TRANSONIC RANGE WITH CONTROLS FIXED 
By Howard F. Matthews 


SUMMARY 

Mach number histories of the motion experienced by a hypothetical, 
small, straight-wing aircraft accelerating at various rates through an 
assumed controls— fixed pitch-down balance change in the transonic range 
were obtained by a differential analyzer. For this particular case, the 
maximum change in normal acceleration is shown to increase with increas- 
ing longitudinal acceleration up to a certain magnitude of longitudinal 
acceleration, after which the maximum change in normal acceleration 
decreases with further increases in longitudinal acceleration. It is 
found that the Mach number variation of the angle of attack for static 
balance determines to a great extent the degree of the Increase in the 
maximum change in normal acceleration. The effects of changes in atti- 
tude, altitude, moment of inertia in pitch, and a factor representative 
of the Mach number range of the balance change on the response of the 
aircraft are also noted. 

Two approximate analytical solutions of the longitudinal equations 
of motion are developed which are based on certain simplifying assump- 
tions indicated by the differential-analyzer results. Examples of each 
of these methods are presented and' the results are shown to compare favor- 
ably with the differential-analyzer solutions. Computation time estimates 
for each method are also given. 


INTRODUCTION 

It is well known that conventional straight-^wing aircraft experience 
unusual balance changes in the supercritical or transonic flight range. 

These balance changes may arise from a loss In lift at constant angle of 
attack, altered wing pitching-moment characteristics, a reduction of 
longitudinal-control effectiveness, or their combined effects. These are 
customarily summarized in reports by a plot of horizontal stabilizer inci- 
dence or elevator angle for steady level flight as a function of Mach 
number . The rapid and sometimes rather large changes in longitudinal con- 
trol deflection necessary for balance shown on such plots, which are 
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usually accompanied by changes in control force for balance, led to the 
belief that a pilot may have difficulty maintaining steady level flight 
as the aircraft accelerates through the transonic range. In addition, he 
possibly may apply adverse control deflection due to his relatively slow 
reaction time and inability to transmit changes in attitude or normal 
acceleration into proper corrective control deflection, especially since 
the effectiveness of the control may be rapidly varying. For these rea- 
sons, the proposal of accelerating through this region with fixed longi- 
tudinal control was advanced, with the further supposition that increasing 
the longitudinal acceleration might reduce considerably the effect of the 
attendant balance changes on the motion of the aircraft. It may be men- 
tioned that the foregoing also has its correlative problem in the accel- 
erating phase of the flight of supersonic guided missiles, including those 
launched from high subsonic speed aircraft, 

A brief study was made, therefore, of the effect of longitudinal 
acceleration on the normal acceleration, particularly the maximum normal 
acceleration, for a hypothetical straight-wing aircraft subjected to an 
assumed pitch-down balance change In the transonic range. Due to the 
complexity imposed by the assumed nonlinear variation of the aerodynamic 
parameters with Mach number, a differential analyzer was used to permit 
rapid simultaneous solution of the three longitudinal equations of motion. 

To make possible the solving of these equations of motion without 
the use of specialized computational equipment, simplified approximate 
analytical procedures were developed based on the results from the differ- 
ential analyzer. These approximate analytical methods also lead to an 
indication of the effect of longitudinal acceleration on the estimation 
of steady-state stability characteristics obtained by the research tech- 
nique of programmed control motions of rocket-powered models. 


A Z 


An 


cl 


NOTATION' 


longitudinal acceleration factor, the ratio of the net aerody- 
namic force in the direction of the relative wind (positive 
when directed forward) to the weight of the airplane ^ari^ 


normal acceleration factor, the ratio of the net aerodynamic 
force perpendicular to the relative wind (positive when 

directed upward) to the weight of the airplane f L — ) 

\ CL trim / ' 


drag coefficient 



lift coefficient 
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airplane lift coefficient for steady level flight ( -2- 

\qs 


r) 


lift— curve slope of complete airplane 

lift-curve slope of the horizontal tail 

pitching-moment coefficient about the center of gravity 
/ pitching moment ^ 




qSc 


V 


pitching-moment coefficient about the center of gravity for 
zero i|., 5 g , and C-^ 


static longitudinal stability parameter 

stabilizer effectiveness parameter 

/Sc_ 

elevator effectiveness parameter ( — — 

Vas e 


(If) 


© 


) 


pitching moment of inertia, sing— feet squared 
Mach number 


(i) 


wing area, square feet 
horizontal— tail area, square feet 
thrust, pounds 
velocity, feet per second 
weight of aircraft, pounds 
velocity of sound, feet per second 
wing span, feet 
wing section chord , feet 


wing mean aerodynamic chord 


fc 


b/ : 
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c dy 


, / 


i, feet 
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gravitational acceleration, feet per second, squared 
pressure altitude, feet 
stabilizer angle 
tail length, feet 
mass of aircraft (&- A, slugs 

. W 

tail efficiency factor 
free-stream dynamic pressure , pounds per square foot 

dynamic pressure at the horizontal tail, pounds per square foot 
time, seconds 

lateral coordinate of wing section chord, feet 

angle of attack 

angle of attack for zero lift 

flight-path angle (8 -a) 

elevator angle 

angle of downwash at the horizontal tail 

ratio of damping of the complete airplane to that of the hori- 
zontal tail 

angle of pitch 

mass density of air, slugs per cubic foot 



All angles are in radians uni ess otherwise noted, A dot (*) or double 
dot (••) above a symbol represents, respectively, the first and second 
derivative with respect to time. The subscripts i and s refer to 
initial and static-balanced conditions, respectively. The term "static 
balanced" is used herein to specify the values of a, C^, or 
which give zero pitching moment in steady (no longitudinal acceleration) 
flight. 
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DH^EEiENTIAL-AKALTZER CAIjCUIATIOHS 
FOR HIPOTHETICAL AIRCRAFT 

Method of Analysis 

Since in the transonic speed range the aerodynamic parameters vhich 
influence the motion of an airplane are generally nonlinear functions of 
Mach number, the simultaneous solution of the longitudinal equations of 
motion appeared too time consuming for a normal step— by— step procedure. 
Therefore, the differential analyzer at the University of California at 
Los Angeles was used to solve these equations of motion. 

A number of solutions were obtained corresponding to different vari- 
ations of the important aerodynamic parameters with Mach number , changes 
in initial altitude and attitude, and several magnitudes of constant 
thrust. This was done with a twofold purpose: first, to show the mechan- 

ics of the motion of an airplane subjected to a balance change in the 
transonic range, particularly with regard to the maximum change in A n 
developed, for different constant thrusts (or, as examination of results 
later shewed, approximately constant longitudinal acceleration A l for 
thrust values of practical interest) ; and, second, to provide information 
relative to a possible simplification of the longitudinal equations of 
motion or other artifice leading to an analytical solution. 

It was realized that in the transonic range the aerodynamic parame- 
ters are not only nonlinear functions of Mach number, but also of angle of 
attack at high lift coefficients. This nonlinearity in angle of attack would 
have offered no difficulty in solving the equations of motion on the differ- 
ential analyzer, but it was neglected partly because it was felt to be of 
secondary importance and partly due to the belief that it could not be 
taken into account in any practical or simple analytical approach. 

For each solution, the airplane was assumed initially to be in 
steady flight at a stated altitude, attitude, and a Mach number of 0.90. 

A given constant thrust was then applied to accelerate the airplane through 
the transonic range. 


Equations of Motion 


If damping in pitch and the lag in downwash at the tail are the only 
nonstatic aerodynamic effects considered, the equations of longitudinal 
motion in wind-axis notation are: 


T cos a — q S C D — W sin 7 = mV 


— T sin a — q S C-r + W cos 7 = -m Y 7 


q Sc C m - 


qs-tW 2 




Y 
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where is assumed equal to C^Cor-a,-)) and. Cm equal to Cm 0 + 

+ ^m§ ^e + Cm,-, Cj,. 

t e T i 

The physical limitations of the differential analyzer made necessary 
the assumption that the mass, moment of inertia, and thrust remained con- 
stant with time. An additional simplification was to consider the cosine 
of the angle of attack equal to unity and the sine equal to the angle in 
radians . 


Assumed Characteristics of Hypothetical Airplane 

The geometric and mass characteristics assumed are: 

S, square feet 150 

S^, square feet . . . 30 

c, feet 5 

Z.J., feet . . . . .". . . . .12.5 

W, pounds 7,500 

Xy, slug— feet squared 10,000 

(unless otherwise noted) 


The variation of aerodynamic parameters with Mach number is shown in 
figure 1. Three variations with Mach number are given in this figure 
for the lift— curve slope and the angle of attack for zero lift. The 
applicable variation is denoted in subsequent figures by the letters a, 
b, or . c. Initial level flight at a pressure altitude of 50,000 feet 
was used unless otherwise indicated. The ratio of the damping of the 
complete airplane to that of the tail £ was assumed to be constant and 
equal to 1 . 15 , and the tail efficiency factor’ n^ equal to 1.0. 

Results and Discussion 

Since the terms "static-balanced angle of attack" and "static-balanced 
normal acceleration" occur frequently in the subsequent discussion, it is 
believed desirable to clarify by two examples the definition given in the 
section Notation. For instance, if tail damping is negligible, the static- 
balanced condition can be thought of as a steady turn of a radius of curva- 
ture just sufficient to develop the static-balanced normal acceleration 
A. ^ and static-balanced angle of attack a g at the particular Mach number 

lander consideration. Likewise, the lift coefficient and angle of attack 
assumed by a free-floating model with the pivot point at the assumed 
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center of gravity are representative of static-balanced conditions. 

Usually the response of an airplane to a disturbance is illustrated 
by a time history of appropriate quantities. For this particular study, 
however, it is believed more meaningful to plot the response as a Mach 
number history, since the basic balance change or forcing function, repre- 
sented by Og or A ns , is Independent of longitudinal acceleration when 
plotted as a function of Mach number. 

As nay be seen from the equations of motion, the instantaneous values 
of the aerodynamic loads and associated moments have “been assumed to he 
completely defined hy the angle of attack of the lifting surface; that is, 
the lag in the build-up of circulation due to unsteady motions has been 
neglected. It follows that any dynamic effect of longitudinal accelera- 
tion must arise through its influence upon the pitching response of the 
aircraft. For the pitch-down balance change defined by characteristics 
(a) of figure 1, figure 2 shows that the dynamic effect of longitudinal 
acceleration occurs as an increasing initial speedwise lag in angle of 
pitch with increasing longitudinal acceleration. 

Iel figure 2 and most of the succeeding figures, the average longi- 
tudinal acceleration corresponding to each value of constant thrust 

is noted. The longitudinal acceleration was found to he nearly constant 
over the Mach number range at high values of thrust, hut Increasing devi- 
ations from a constant magnitude occurred at lover thrusts. It is partly 
for this reason that data at low thrusts are omitted and partly due to the 
fact that the extreme changes in attitude and altitude and the relatively 
large time interval necessary to traverse the balance change associated 
with low thrusts preclude their practicability. 

In figure 3 are plotted Mach number histories of the angle of 
attack and normal acceleration for the same average longitudinal acceler- 
ations as in figure 2 and for the static-balanced condition. It Is appar- 
ent from this figure that at a relatively low value of the response 

of the aircraft follows closely the static-balanced condition with increas- 
ing departures occurring at higher accelerations. At some magnitude of 
A? , a peak in the maximum change of normal acceleration is reached, 

which, for the conditions assumed, is appreciably greater percentagewise 

than for a low value' of , although the actual nagnitude of the 

a v 

increase in g Is small due to the extreme altitude. This effect of 
longitudinal acceleration upon the maximum change in normal acceleration 
(or to the ratio of the naximum change in dynamic to static-balanced nor- 
mal acceleration, which, for convenience, will he referred to hereafter 
as the normal-acceleration response ratio) is largely dependent upon the 
variation of the static— balanced angle of attack with Mach number. 

Figure k shows the results of a number of solutions hy the differential 
analyzer for three variations of a g but an identical variation of A^ 

with Mach number. The variations in a g were obtained hy changing the 
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Mach number variation of a Q and C-^. It is seen that the most pro- 
nounced effect of longitudinal acceleration upon normal acceleration 
occurs for the condition of the greatest variation of <x 8 with Mach 
number . 

Comparatively large changes in initial flight-path angle (0° and 30°) , 
moment of inertia (5,000, 10,000, and 20,000 slug— ft 2 ), and altitude 
(30,000, 40 , 000 , and 50,000 ft) have relatively little effect on the mag- 
nitude of the peak of the normal-acceleration response ratio. A change 
in Mach number width of the balance change, the width being defined herein 
as the Mach number increment between the initial Mach number and that for 
maximum change in A_, , also has little effect on the peak normal— 

acceleration response ratio if A^ and a g plotted as a function of 

percent of the Mach number width of the balance change remain the same. 

As illustrated in figure 5 , independent changes in moment of inertia, 

Mach number balance change width, and altitude have an appreciable effect 
on the magnitude of the average longitudinal acceleration at which the 
peak normal-acceleration response ratio occurs. The assumed balance 
change in this instance is of a larger magnitude than that for previous 
figures, but it is believed that the results in figure 5 indicate the 
general trend. 

From the data presented, it is apparent that, if the Mach number 
width of the balance change is of the order of 0.05, the total thrust 
capabilities of transonic research aircraft may place them in a region 
where the normal acceleration is appreciably affected by longitudinal 
acceleration. 


SIMPLIFIED COMPUTATION METHODS 

Since a differential analyzer is not always available, the develop- 
ment of relatively quick approximate analytical methods by which the 
response of an aircraft to longitudinal acceleration could be computed 
was believed desirable. In the following two methods, the validity of 
the assumptions made in developing these methods was indicated by the 
results from the differential analyzer and was substantiated in each case 
by comparing the normal acceleration computed by the analytical method 
with the differential-analyzer answers for two values of longitudinal 
acceleration. These methods are based on reducing the three longitudinal 
equations of motion to a single differential equation in angle of attack 
which is of the second order, linear, and nonhomogeneous . 


Modified Step by Step 

Theory.— Examination of Mach number histories for the longitudinal 
accelerations considered revealed that the flight-path angle (but not 7) 
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was relatively constant. Since the time interval to traverse the "balance 
change is small, resulting in little change in altitude, the density of 
the air and speed of sound also can he considered unchanged during the 
flight. 

If V is replaced "by aM and V hy sM, the two force equations 
may "be combined by eliminating thrust and the longitudinal equations 
reduce to the following: 

“ “ 4 - [ ( 2 * <0ltt + ° D) + (* + f sln n) J ] “ + 

K 

( a ^) K M ^ + (f coe 7i X K (1) 


& ~ ( 2 M 2 1" Cjjiq + it + C m g e 8 e +■ C ma (ar-a 0 ) 1 — 

• T K L J 

( 1 ^) («*£*) < 2 > 

K 

where the terms enclosed with parenthesis with a subscript K are con- 
stant. Differentiation of equation (l) gives 


a = 

»-[d 

r *?r) M < c i« + c »> + 

j'* + S. sin 

71 

) h] 

4 - 



■K 

\ 


K J 



( 1 paS 
_ \2 m , 

) MCc^ + c D ) + (vl 

+ — sin 7 i 
a 

) 

a): 

a + 



K 


K 



/ 1 paS 

V 2 m 

) M a o + ({ cos 

r 0 (*) 






K 

K 




It is advantageous to evaluate tiie 

term(l) 

at 

this time. 
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dt M 2 


Substituting equations (l), (2), and (4) into equation (3) to eliminate 
0 and its derivatives and noting that 


— Cm a. = C_ + C„ i+ + Cm 8_ — . Cm 

“a s m o “4^. t m 6 e ® “a ° 

the equations of motion reduce to the following single, second— order, 
linear, nonhomogeneous differential equation with variable coefficients: 


'd + ba + ka = f 


(5) 


where 


15 - (i tt) m + c d) -(« + § “1“ ?i) H + 

K K 




M nt C 


^t 


0 + £) 


“ <Cl - + Cl,) 


0*-s-O i]*G¥*) 

K K 

( M 2 + | M sin 


CD)] 


M 2 
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= _ (l oa. s Sc \ 

V2 i T yfc 


: ™o, “s + 


1 ^ 


L 1 M nt 5 [(§ 


MC Itt a b + 




Equation (5) is recognizable as the fundamental equation of simple, damped, 
vibratory motion with, thq exception that in this instance the coefficients 
b and k, representative of the damping and the spring constant, respec- 
tively, are nonlinear functions of Mach, number or time. Therefore, consider 
the case in which the Mach number range under investigation is divided into a 
number of intervals, the interval being chosen sufficiently small .so that the 
variables b and k can be assumed constant in magnitude and f a linear 
function of Mach number. The solution of equation (5) for constant coeffi- 
cients can take any one of three forms, depending on whether both roots of 
the auxiliary equation are real and unequal (b/2) 2 >k, real and equal 
'(b/2) 2 = k, or complex k> (b/ 2) 2 and can be found in any standard text 
on differential equations or vibrations. The following is the solution 
(in terms of M since t = for the most common case of k>(b/2) 2 : 
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where 

f Is represented by f _ + — AM 

AMO dM 



F 


0 


f AMO ~ 
k 


A = a AM=0 - F 0 


^AMO + — ~ 

*/k - (b/2) 2 


The Mach number interval is denoted by AM and the initial 0^=0 
is obtained from equation (l) by assuming Q equal to zero. The Mach num- 
ber history of angle of attack of the aircraft is computed step by step by 
means of equations (6) and (7) for each Mach number interval, the initial 
conditions of one interval being equal to the end boundary conditions of the 
previous interval. The motion of the airplane can then be completely defined 
by obtaining the angle of pitch from equation (l). 

Results and discussion.— The necessary curve fitting and computations 
of. an example are shown in figure 6 and table I, respectively. The selec- 
tion of AM (a constant value of 0.02 for the example shown) is obviously 
an important factor and is a compromise between the AM for minimum com- 
putation time and the AM for greatest accuracy. The optimum value is 
clearly that magnitude, which will give the desired accuracy in the minimum 
time and may be best selected by the examination of the computations at a 
number of values of AM. The AM need not be the same for each interval, 
but the computations are reduced slightly if it is constant. The results 
of this example and .one at a higher value of A are compared with the 

4 a,v 

differential-analyzer answers in figure J. It is evident that by this 
method good correlation is obtained with the differential— analyzer Mach 
number histories. 

Initially, it may appear that this method is at least as time 
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consuming as the normal step-by— step procedure, hut a careful scrutiny of 
the equations defining bj k, and f shows that only portions of these 
factors are affected by M, and then only in a simple manner. This is a 
very important time-saving characteristic if Mach number histories are 
desired for a number of longitudinal accelerations. Some simplification 
in the factors b, k, and f usually may be made. For instance, Cp 
is normally negligible compared to even in the transonic region, 

and the portion of k containing the tail damping may be negligible, 
depending upon the amount of static stability. It has been estimated 
that, given the necessary aerodynamic data in the form as shown in figure 
1, a Mach number history in angle of attack for one longitudinal accelera- 
tion would require approximately two and a half computer days and that 
each succeeding longitudinal acceleration would necessitate about one 
additional computer day. 


A study of the terms forming the coefficients b, k, and f of 
equation (5) indicates that there are two basically different effects 
of longitudinal acceleration; one that arises from the terms which are 
functions of ft; and the other from the fact that different longitudinal 
accelerations cause the assumed disturbance to be traversed in different 
lengths of time or periods. The latter effect is analogous to that of 
first— order vibrating system in which the amplitude of the response 
depends both on the magnitude of the input disturbance and on the ratio 
of its frequency to the natural frequency of the system. For simple 
disturbances then, an estimation of the longitudinal acceleration for 
maximum change in normal acceleration may be made by computing the ft 
which will result in the period of the disturbance being equal to the 
average natural period of the airplane. 


A method of research in the transonic range which is being increas- 
ingly utilized 1s the rocket— powered model flight, wherein the controls 
are programmed to have simple motions throughout the flight and the 
response of the model measured by appropriate ins truments - Much useful 
information is obtained from these tests, such as an estimation of the 
static and dynamic— stability characteristics derived from the period and 
damping of the oscillation following the control deflection. The period 


\ fc - (t/2)' 


and time to damp are inversely proportional to the coefficients 
and b, respectively, of equation (5). It appears from an examination of 
these coefficients that appreciable errors can occur in the determination 
of the stability of the full-size configuration from the powered portion 
of such flights of high longitudinal accelerations are developed during 
this period and for this reason the coasting portion of flight is used for 
the data analysis. 


E qulvalent— Di s turbanc e Method 

Theory .— The differential-analyzer results indicated that the Mach 
number history of (or Ct ) and <x s had considerable effect on the 

response of the aircraft during longitudinally accelerated flight. This 
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fact suggested the concept of an "equivalent disturbance"; that is, one 
in vhlch the Mach number variation of the aerodynamic parameters which 
define C Tg and a g could be changed to suit same particular purpose, 

and, provided the Mach number hiBtory of C L and a g re m ai n ed the 
same, the response of the airplane would remain relatively unaltered. 


A close approximation to the correct solution is then obtainable by 
assuming C^, C^, n^. Cj^, de/da, £, and M (where it appears in 

the equations of motion) constant over the Mach number range and adjusting 
it + Cn R 8 e and a Q so that the Mach number variation of 


C wi Q + Cm^ 

3, and a_ 


remains the same. 


For this method, it is convenient to 


rewrite the equations of motion in terms of the angle of attack, the angle 
of zero lift, and the static angle of-attack, all measured from initial 
condition. With C D neglected, equation (l) can then be rewritten as 


a 




m 




( 1 £ 5 §. 

\2 m 



+ 


■o) 



( 8 ) 


-where, as before, the terms enclosed within parenthesis or brackets with 
a subscript K are constant. For steady— state conditions at t=0, 
equation (8) becomes 


K K 

( 8 cos 7i\ 

v a M 4 


a °t=o + 


(9) 


Subtracting equation ( 9 ) from (8) gives 
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(I), 


M , 

M J ^=° 
E 


( 10 ) 


Similarly, the p itching-moment equation becomes 


(— .)-(l^H0 s ) [«, K e 

K K 

(=) 4 


(id 


Combining equations (10) and (11) as before, the resulting differential 
equation is 


Ax + b f Ax + k f Ax = f 1 


( 12 ) 


vhere 


b* = 


1 MS i, c , * + f Sia , 1 paS^t 2 ** 

[2^ MC Iu M“ 2 ~~1Z 


c I«ot 


(*s)] 


E 


k* - f"- i — — M S Cjj, + i paStZt - M ixt C L g (i M c 

J_ 2 I y % 2 i T ^ b \Q m 


M + | Bln y ± 
M 


)] 


E 


f » 
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K 


the negligible constant term 


1 P^t^t 2 . \ 

2 ~ y nt C S; ^ Ma ^=0j K 


"being omitted from f*. Equation (12) is now of the form where the prin- 
ciple of superposition applies and Duhamel t s integral theorem may be used 
to calculate Zsx due to the irregular variation of the disturbance f*. 

A semigraphical application of this theorem which was applied in reference 
1 to the calculation of the motion of an airplane under the influence of 
irregular disturbances and more recently in reference 2 was -used for the 
present computations . This method reduces essentially to the solution of 
equation (12) for a unit step— forcing function f ' (l) and the combination 
of this response with the variable f * by means of Duhamel’s integral 
through use of a simple graphical procedure. 


Results and discussion.— The assumed equivalent-disturbance param- 
eters , along with the original variation of the parameters with Mach num- 
ber for the same examples as shown in figure 7 are presented in figure 8. 
The values of C^, and de/da represent average 

values over the range from the initial Mach number to the Mach number for 
maximum change in A^; it was found that these magnitudes give good 


results. The Mach number' at the maximum change in A n D was used in com- 
puting the coefficients b*, k*, and f*. 


Figure 9 illustrates the graphical integration procedure for one of 
the examples at a specified instant of time, and is discussed in detail 
in the appendix. 


Figure 10 shows the good agreement between the computations by this 
method and the differential-analyzer answers for the same two longitudinal 
accelerations used in the modified step-by— step method comparison. The 
method is relatively rapid, the first case requiring approximately two 
computer days and each succeeding case about three-quarters of a computer 
day. It should be noted that if- the effect of M upon b’ and k* is 
small, the response to a step input f’(l) need be computed only once, a 
very important time-saving factor if a number of longitudinal accelerations 
are to be investigated. 


An evident inaccuracy arises from the use of a constant Mach number 
in the simplified equations of motion in place of the true Mach number 
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variation over the range "being investigated. It Is believed, therefore, 
that this method should be applied only to a balance change, the Mach 
number range of which Is small. 


CONCLUSIONS 

As a result of an investigation to determine the effect on the 
motion experienced by a hypothetical, small, straight— wing aircraft when 
accelerated at various rates through an assumed controls— fixed pitch- 
down balance change in the transonic range, the following conclusions 
may be made: 

1. The maximum change In normal acceleration increases with increas- 
ing longitudinal acceleration up to a certain magnitude of longitudinal 
acceleration, after which it decreases with further increases In longi- 
tudinal acceleration. 

2. The Mach number variation of the angle of attack for static bal- 
ance determines to a great extent the degree of the effect of longitudinal 
acceleration on the normal acceleration. 

3. Changes in altitude, initial flight-path angle, moment of inertia 
and a factor representative of the Mach number range of the balance change 
(defined herein as the Mach number width) have little effect on the magni- 
tude of the peak of the ratio of the maximum change in dynamic to static- 
balanced normal acceleration, but an increase in moment of inertia and 
altitude and a decrease in Mach number width decreases the magnitude of the 
average longitudinal acceleration at which this peak occurs. 

k. The two approximate analytical methods for solving the equations 
of motion which were developed give good results for the assumed balance 
cha n ge. The modified step-by-step method appears to be applicable to all 
problems of this type. However, the accuracy of the equivalent-disturbance 
method has not been adequately Investigated so Its use should be limited to 
cases similar to the one Investigated. 

5. The methods presented herein may be used to study the effects 
of longitudinal acceleration on the results obtained from rocket models 
flown to determine the static and dynamic stability characteristics and 
as an aid in planning flight programs to avoid the effects of longitudinal 
acceleration on these characteristics. 

It is difficult to generalize on the basis of the results of this 
brief investigation, but It is believed that the above noted trends 
indicated should apply to any similar balance change. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif. 
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APPENDIX 

CALCULATION OP THE INCREMENTAL ANGLE OF ATTACK Ax 
FOR ANT VARIATION OF TEE FORCING FUNCTION f ' 

The solutions of equation (12) for a unit step function f*(l) can 
take any one of three forms, as explained in the section on the modified 
step— by-etep method, hut only the solution for k f >(by2) 2 will he pre- 
sented in the following equation* 


Axk 1 

f*(D 


1 + e -(Wa) t ^ , -(t«/2) 

J k*-(h*/2) 2 


cos 


J k*-(b'/2) 2 t 1 


sin 



The graphical method for combining the response of a unit disturb- 
ance with the variable forcing function f' can best be understood by 
following a sample computation. The work sheet for an example is shown 
in figure 9, and the procedure is as follows: 

1. Plot Axk*/f ’ (1) as a function of time to some convenient 
scale. 

2. Plot f*. as a function of- time using the same time scale. 

3. Select a time t Q at which the Ax due to f* is desired. 
Project the point on the f * curve corresponding to this time horizon- 
tally until It intersects the 45° line. This line is then deflected 
vertically until it intersects the horizontal projection of Axk , /f‘ , (^) 
at t=0. This establishes the point labeled (1). 

4. The ordinate of the f * . curve at same time less than t Q by 
the amount At is next projected as before until It intersects the 
horizontal projection of-the Axk'/f'tl) at At. The point (2) is 
obtained in this manner. 

5. Other points are similarly obtained to complete the closed 
curve for the time t D . Note that the addition of t on the f’ curve 
and t on the Axk*/f*(l) always equals t D . The area encompassed by 
the curve is proportional to Ax at time t Q . 

6. The area is found by integrating in the direction (0), (l), 

(2), etc. If a counterclockwise path is followed in enclosing the 
area, the value is positive regardless of the quadrants involved and 
vice versa for clockwise integration. 
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7. Other curves are drawn for different times. The final step is to 
correct the areas into the associated with the function f * "by multi- 

plying the areas by the appropriate scale factors of k*/f r (l) and f* 
and dividing by k* . The lift coefficient is equal to Ct_ 

"^assumed 


(c^-n + Ax — a_ ) . The correct angle of attack can be obtained from 


u assumed 

Ct 

the formula a = 1 


and the angle of pitch from equation Cl) 


11 true 
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TABLE I.- COMPUTATIONS OF NORMAL ACCELERATION BY MODIFIED STEP-BY-STEP METHOD 
[Characteristics (a) of figure 1, average longitudinal acceleration of 1.373 
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Figure /.- Assumed variation of aerodynamic parameters with Mach 
number for the hypothetical airplane . 
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Mach number, M 


Figure 2. - Effect of longitudinal acceleration on the initial 
response in pitch. 



Normal acceleration factor , A„ Angle of attack 
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Figure 3.- Effect of longitudinal acceleration on the angle of 
attack and normal acceleration factor . 
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Figure 4- Effect of longitudinal acceleration on the normal 

acceleration response ratio for three variations of a s with 

Mach number but identical variation of A„ with Mach 

"s 

number. 
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Average longitudinal acceleration, A tgv , at the 
peak normal acceleration response ratio 

Figure 5.- Effect of changes in moment of inertia, Mach number 
disturbance width, and pressure altitude on the average 
longitudinal acceleration at which the j>eak normal acceler- 
ation response ratio occurs. 
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Figure 6.- Example bf curve- fitting method. Characteristics 
(a) of figure f, average longitudinal acceleration of t.37 . 











y 

s 








s 

— 

* 








s' 





— - 



■* “ 






- 


1 





1 1 






NACA RM A9J2 6 


27 


Modified 

step-by-step Differential 
method analyzer Thrust, its Aj 

12,500 f.37 

40,000 5.00 

Characteristics (a) of figure i 



.90 .94 .98 i.02 1.06 UO 

Mach number, M 


Figure 7. - Comparison of normal acceleration factor computed 
by the modified step-by-step method and by the 
differential analyzer. 
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Mach number , M 


Figure 8.- Mach number histories of original and equivalent- disturbance 
method data. Average longitudinal acceleration of L37 . 
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Figure 9. - Example of the graphical procedure for the determination of the incremental angle 
of attack A a. Average longitudinal acceleration of 1.37 . 
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disturbance Differential 
method analyzer Thrust , lbs Aj 

12,500 ! .37 

40,000 5.00 

Characteristics (a) of figure / 
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Mach number, M 

Figure /O - Comparison of normal acceleration factor computed 
by the equivalent- disturbance method and by the differ- 
entia! analyzer. 
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